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Abstract 

Assumed that the parameters of a generahzed hypergeometric function depend hn- 
early on a small variable e, the successive derivatives of the function with respect 
to that small variable are evaluated at e = to obtain the coefficients of the e- 
expansion of the function. The procedure, quite naive, benefits from simple explicit 
expressions of the derivatives, to any order, of the Pochhammer and reciprocal 
Pochhammer symbols with respect to their argument. The algorithm may be used 
algebraically, irrespective of the values of the parameters. It reproduces the exact 
results obtained by other authors in cases of especially simple parameters. Imple- 
mented numerically, the procedure improves considerably the numerical expansions 
given by other methods. 
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1 Introduction 



Special functions of tlie hypergeometric class have deserved a considerable 
interest in the last decade due to their connection with Feynman integrals in 
quantum field theory. Recurrence relations among hypergeometric functions 
allow one to reduce the number of Feynman integrals to be computed in a given 
process. In reciprocity, calculations with Feynman integrals have revealed new 
relations between hypergeometric functions pLj. 

Dimensional regularization of Feynman integrals can be conveniently carried 
out if one knows the e-expansion of the related generalized hypergeometric 
functions, that is, an expansion in powers of a small variable e on which the 
parameters of the function depend linearly. Several procedures have been sug- 
gested to deal with this issue. Algebraic methods applicable in the cases of 
integer, half integer or rational parameters [2|l3|H|5|6|rf|l8|l9] lead to expansions 
whose coefficients can be written in terms of generalized polylogarithms. Com- 
puter packages implementing algorithms based on the Hopf algebra of nested 
sums [rUfTT] are available. Let us mention, for instance, a C-I--I- program [12] 
in the framework of GiNaC, XSummer [T3], which uses the computer alge- 
bra system Form, HypExp and HypExp 2 [Tip5] . based on Mathematica, and 
HYPERDIRE |16], useful for differential reduction of hypergeometric functions. 
A recently proposed package, NumExp [18], with a different strategy, allows 
one to get the e-expansion as a Laurent series whose coefficients are evaluated 
numerically by using a multi-precision finite difference method. 

Here we are concerned with the e-expansion, to all orders, of generalized hy- 
pergeometric functions of one variable, pFg, or of several variables, like the 
Appell and the Kampe de Feriet functions [T7], with parameters depending 
linearly on e. No restriction is placed on the parameters, that may be complex. 
Our procedure stems from the well known series expansions of those functions 
and, similarly to the method in Ref . [18] , consists in obtaining the coefficients 
of the Laurent series by multiple derivation with respect to e followed of par- 
ticularization for e = 0. Instead of using (approximate) numerical derivation, 
like in Ref. [IB], we benefit from simple formulae for the derivatives, of any 
order, of the Pochhammer and reciprocal Pochhammer symbols, which allow 
us to obtain explicit algebraic expressions for the successive terms of the e- 
expansion. The validity of our expressions is subordinate to the convergence of 
the series expansions of the functions considered. We assume in what follows 
that the values of the variables are such that the needed convergence is guar- 
anteed. Otherwise, such expressions have a purely formal character. Needless 
to say, analytic continuation allows to cover a wider range of values of the 
variables. 

We present in Section 2 our algorithm. For the sake of clarity, we consider the 
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£-expansion of 



QJl, . . . , Oip 




Z 



(1) 



with upper and lower parameters of the form 



ai = Ai + aie, i = l,...,p, 



l3j = Bj + bje, j = l,...,g, 



(2) 



without any restriction on the values of Ai, ai, Bj, bj. The procedure can 
be trivially extended to Appell and Kampe de Feriet functions. By way of 
illustration, wc report in Section 3 the results of applying our method to sev- 
eral examples of generalized hypergeometric functions found in the literature. 
Section 4 contains some pertinent comments. We have devoted Appendices A 
and B to provide with expressions for the derivatives of the Pochhammer and 
reciprocal Pochhammer symbols. 



2 The algorithm 

As said before, our starting point is the series expansion 



where one uses the Pochhammer symbols, = r{a + n)/T{a). 

We say that all lower parameters are regular if all Bj {j — 1, . . . ,q) are 
different from zero or a negative integer. In this case the £-expansion is a 
Taylor one, free of negative powers of e. If one or several lower parameters 
are singular, that is, if some of the Bj are zero or a negative integer, negative 
powers of e may be present and the expansion is a Laurent one. We treat in 
turn the two cases. 




(3) 
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2.1 All lower parameters are regular 



Let us denote by Cn{z) the coefficients of the e-expansion 
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Obviously, 
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Let us introduce the abbreviations 



k\ da'' 



1 1 



(6) 



to represent the derivatives of the Pochhammer and the reciprocal Pochhamer 
symbols with respect to their arguments. Then, derivation with respect to e 
in Eq. ([3]) gives 



& 



de 



n P^l 



Oil, ■ ■ ■ , Oip 

^/3l,...,/3g 



oo n 



z =n 



ni 



m=0 A:i=0 



n„_i 



X 



k2=0 kp=0 
rip "p+i 



A;p+i=0 



kp+q-l=0 

where we have denoted 
ni = n — ki — k2 — ■ ■ 



ki, (/ = l,2,...,p + g-l). 



(7) 



The expression of C„(z) then results immediately from the right hand side 
of Eq. ([7]), by suppressing the initial factor n! and substituting Ai for ai 
{i = 1,2, ... ,p) and Bj for 13 j {j = 1,2, ... ,q). Simple explicit expressions 
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for V^\Ai) and Q(^^(i?j) are obtained in Appendices A and B, respectively. 
They are 



V^''\A,) = 6k,o, V^\A,) = for k>m>0, 



(9) 



m—k 



V^J:\A,) = i-ir-'^ Y: (-1)' T k) iA,)i for m>k, (10) 



1=0 
(0), 



m— 1 



/!(m-l-0! {Bj + l)>'+^ 



(11) 

for m, A; > 0, (12) 



The s{n, k) in Eq. ffTOl) represent the Stirhng numbers of the first kind [HI 
§26.8], whose generating relation is 



[ln(l + x)]'' = A;! V s(n,fc):^, 



\x\ < 1. 



^.^ Some lower parameters are singular 



Let us consider now the case of . . ., /3r (0 < r < g) being singular because, 
for j = 1, . . . , r. 



S,- = —Nj, Ni < . . . < Nr nonnegative integers. 



(13) 



The series expansion in Eq. can be written in the form 



1 , . . . , Oip 

. . . ,l3g 



^'^^0 (/3iV---(/3,)™^ 



+ 



-1 N2 



E 



+ 



E 



pjm 



+ . . . 



+ 



E 



(ai)r 



p )m 



blb2 ■■■br „=jv. + l (/3l)^ ■ ■ ■ mjPr+l)m " " " (/3g)n^ ^' 



(14) 



where a sum is void if the lower limit of the summation index is larger than 
the upper one. We have represented by l/(/3j)^ the "regularized" reciprocal 
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Pochhammer symbol 
1 _ bjE 



p. = -Nj + bj e, < Nj < m, 



(15) 



whose derivatives with respect to its variable, 



(16) 



are calculated in Appendix B. The result is 

(-1)' 



m—l 
/=0, lytNj 



Z!(m-1 - /)! (/3j + 0*^+1 



(17) 



and, for e = 0, 



m—l 



e=0 



The e-expansion is now of the form 





• • ; Q^p 

















E 



(19) 



The coefficients C„(z) are immediately obtained as a sum of the corresponding 
coefficients of the e-expansions of the sums in the right hand side of Eq. (|T4l) . 
which can be evaluated by means of the algorithm used to obtain Eq. ([7]) from 
Eq. 0. 



3 Some examples 



The procedure sketched above may be used to get algebraic expansions, when 
the parameters Ai and Bj take especially simple values. Nevertheless, we find 
our method most suited to obtain numerical expansions, for arbitrary values 
of the parameters, by means of an implementation of the procedure in FOR- 
TRAN or Mathematica or any similar language. The algorithm is so simple 
that every user can construct an efficient package. In this section we consider 
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several examples for which results, obtained by different methods, are avail- 
able in the literature. Our procedure is able to reproduce the exact algebraic 
expansions and improves considerably the numerical ones. 



3.1 Algebraic expansion of a Gauss hypergeometric function 



As an example of the algebraic use of our algorithm, let us try to obtain the 
first terms of the ^-expansion of the Gauss hypergeometric function 



2F1 



di 6, Q.2 S 



(20) 



According to Eq. ([7]) and its form for e = 0, the coefficients of the e-expansion 
are given by 



00 m n 



n—ki 



m=0 fci=0 



fc2=0 



In the particular case of Eq. (!20|) . one has Ai = 0, = 0, Bi 
values of the parameters, Eqs. ([9]) to (fT2|) reduce to 



(21) 

1. For these 



^^(0) = ^^, ^i°H0) = 5„^,o, Pi'H0)=0 for fc>m, (22) 



pW(o) = (_l)™-^-s(rn,fc) for m>A;>0, 



(23) 



1x^1 Hr, 



m\ TT[l ml 



Q('=)ri) = (-1)^= y ^ ' - 



for m, /c > , 



where one uses the harmonic number if^, 

i^m = E 7 • 

/ 1 



(25) 



(26) 



Then, the first coefficients of the ^-expansion are 



00 



Co(^) = E^(^i°Ho)) Qt\i) 



m=0 



(27) 
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m=0 

0, 



(28) 



oo 



+ {a,a,{v^\0)y + {al + al)V^^\0)V^\0)^ Q^\l) 
= ^ — aia2 (s(m, 1)) 



m=0 



m! 



m! 



which, by using the representation of the Stirhng number 

s(m,l) = (-1)™-^ (m-1)!, 



(29) 



(30) 



and the definition of the dilogarithm function, Li2(^) [T9| 25.12.1], can be 
written in the form 



C2{z) = ai a2 Li2(z). 



(31) 



oo r 



+ &1 («i «2 [V^r^^ (0)) ' + (a? + al) P W (0) (0)) Q^;^) (1) 
+ ((a?a2 + aia^)P«(0)Pi^)(0) + {a\ + 4) Vi^\0) r^\0)) Q^^\l) 

oo m 2 

E (P£H0)) Q«(l) 

oo m 

+ E -Tai«2(ai + a2)P£H0)Pi^)(0)QW(l) 



m=l 



m=2 



oo 



m=l 



Hr, 



ml 



oo 



-(01 + 02) E — r -^l"^' 1) -^l"^' 2) — : 

^2 ^■ 



which, by using the representations Eq. ( I30l) and 
s(m,2) = (-l)™(m-l)!if„_i 



(32) 



(33) 



of the Stirhng numbers, can be written in the form 
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The last equation can be expressed in terms of the polylogarithm function, 
Lin(2;) [ini Eq. 25.12.10], and the Nielsen's generalized polylogarithms, Sn,p{z) 
[20], to give 

Cs^z) = ai a2 [-bi 113(2;) + (ai + 02 - bi)Si^2{z))]. (35) 

We see that our method reproduces the results obtained by Kalmykov [5', Eq. 
4.7] by help of relations presented in Ref. [21]. 



3.2 Algebraic expansion of an Appell function 



Del Duca et al. [22^ Eq. (5.19)] have considered the algebraic e-expansion of 
the Appell function 



1, 1 + e 
l + e,l + e 



Xi, X2 



(36) 



by using an algorithm based on the algebra properties of nested harmonic 
sums. From the series expansion of the Appell function 



Xi, X2 



00 00 

E E 

mi=0 m,2=0 



7712 

2 



{Pl)m, (A 



2)m2 



mi! 7712! 



(37) 



application of our method would give for the coefficient of in the e-expansion 
the expression 



00 ao m\ m2 n n—ki 

Cn{xi,x2)= E E E E a'2^p!^:LM2) 

mi=0 m2=0 "^1- ■"^2- ^^=9 A;2=0 
n— fci— ^2 

X E Qtf{Bi) b^-^'-''^-^^ Qt;''~''~'''\B2),i38) 

fc3=0 

which, for the values of the parameters in Eq. (l36l) reduces to 
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mi=0 m2=0 "^1- "^2- fc^^g 
n— A;2 

X E Q^^'J^l) Qf^-'^^^l), (39) 

fc3=0 

By using Eqs. (ED and ([lO]) with Ai = 1 and Eqs. ([24]) and (ESD, we obtain for 
the first coefficients of the e-expansion (understanding that a sum is zero if 
the lower hmit of the summation index is larger than the upper one) 



nil 1^2 
mi=0 m2=0 "^1- "^2! 



mi m2 



mi=0m2=0 \ J 

■F4l,l-l,l-xi,X2). (40) 



00 00 



Ci{xi,x2)= J2 E 



™mi m2 
Xi X2 



)(0) 



p(0) ^^ 
' mi+m2 K-^J 



(1) 

mi+m2 ' 



(1) qS(i) q^^Jw 



^iiVm2(i) + 



CXD 00 

E E 

mi=0 m2=0 



nil + ^^2 

, mi , 



--mi 
2 



•^1 ^ "^2 ^ [-^mi+m2 -^^mi H^n^ ■ (41) 



oo oo 



C2(Xi,X2)= E E 



mi m2 
^1 ^2 -p(O) 



mi+m2 
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^s+m2(i) 



+ <^m2(i)(QS(i) 2^^1(1) + 
+ <^m2(i)(QS(l) + 2S(l) 2^(1) + 



00 00 



nil + ^^2 



mi=0m2=0 V 



(-1) 



mi+m2 



mi m2 
1 2 



mi 



mi+m2 



E 



Hmi+m2 (-f^mi + ^^2) + I ^^1(^1) 



Z=2 

1 



m.2 



E(-i) 



(42) 



The same coefficients, with a different notation, have been obtained by Del 
Duca et al. [221 Eq. (5.19)]. The reduction of the double sums (in mi and 
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1712) to known functions does not seem to be an easy task. In \22\ Eq. (5.90)] 
transcendental functions Ai{i,j,k]Xi,X2) have been introduced to represent 
the different terms of the coefficients of the e-expansion. The properties of 
those M. functions are discussed in an Appendix of the same paper. 



3.3 Numerical expansion of several generalized hypergeometric functions 



As a first example of numerical application of our algorithm, we consider the 
expansion of 



4^3 



(43) 



whose first terms have been obtained by Carter and Heinrich [23j by using 
the HypExp package [T^fTS] and by Huang and Liu with their own NumExp 
package [18] . Our Table 1 shows their results together with those obtained by 
means of a pedestrian Mathematica implementation of our method. Obviously, 
the m-sum in the series expansion of the hypergeometric function must be 
truncated at, say, m = M. The value of M determines the accuracy of our 
results, that may be controlled by successive increments of M. Notice that, 
as all parameters are rational, Mathematica gives, for the truncated m-sum, 
exact results in the form of rational numbers. The values shown in Table 1 are 
the decimal approximation of those rational results. For the parameters and 
variable in Eq. (143|) a value of M = 50 is enough to give the digits reported. 



In the second example, 



5^4 



^e, —e, —3e, —5e, —7e 
^e, 45, Qe, 8e 



(44) 



we realize that all lower parameters are singular, as 
Bi = B2 = B3 = B4 = . 



(45) 



Equation (fT4|) adopts then the form 



5^4 



''ai, . . . , as 
\Pi, . . . ,/34 



&1&2&3&4 



m=l (/3l) 



E 



(46) 
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Table 1 

Coefficients of the e-expansion of the generahzed hyper geometric function given in 
Eg. iSD. 



e order 


HypExp 


NumExp 


our method 


e° 


1 


0.999999999999999999 


1 




-4.27969 


-4.27968776167885 


-4.27968776167886 




-26.6976 


-26.6975474079 


-26.6975474079466 






195.87119 


195.871193504205 






-7313.7 


-7313.74176765086 








90693.2356441548 








-1426862.01660383 








17612046.1433229 








-233969019.148423 


£9 






2846673719.75988 








-35635855655.1898 



Nevertheless, this particular case is especially simple, as 

Ai = A2 = A^ = Ai = A^ = 0, (47) 



and the negative exponent terms of the Laurent expansion become identically 
equal to zero. Equivalent to the expansion in Eq. fH6|) . we may write the 
following one 



5^4 



'^ttl, . . . , ^5 

. . . , /34 



1 + 



(2]^ (24(25 6. 



E 

m=l 



m— 1 



;i+/3 



lim-l 



(l + /34)m~i rn\ 



which can be treated as described in the case of no singular lower parameters. 
The algebraic procedure gives for the first coefficients of the e-expansion of 



5^4 



(X\£^ • • • ; Oi^E 
biE, . . . , 645 



the expressions 



Co(z) = l, 
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Table 2 

Coefficients of the e-expansion of the generahzed hyper geometric function given in 
Eg. (iMD. 



P OTTlPr 


T~rvnl^^i"V"n IVn tti T^^i'vn 

J. X y t ^ 1 j ^ L/ J. > 1 4 1 1 1 1 J y\ L/ 


miv TTipthnn 


e° 

C 


1 n qqQqqqQqqqQqQQqqqqq 

X \J , U U u U u U u U u U u U U U U ij \J \J ij 




C 




1 89,5324321 843fi 




—2 299043 —2 2990427423 


—2 29904274238202 


C 


4fi902 55 4fi901 9 


5,5 4fi901 90380554 




— 1014 39 —1014 39 


— 1014 392422652,34 


C 




1 5729 382951 742 






-216907.17756543 






2728106.3284185 






-31818216.953372 


£9 




348410894.51286 


£lO 




-3617363078.8137 


Ci(z) 


01020.30.405 / , XX 

= , , , , ( ln(l 2)) , 
6162O3O4 






= (01 + . . . + 05 bi ... 64) „ 
Ol • • ■ O4 z 


[ln(l - z)f , 



and so on. The numerical results of our method, for the values of the param- 
eters and the variable given in Eq. (jH]), are shown in Table 2, together with 
those obtained by using HypExp [23] and NumExp [_18j. 

Huang and Liu [H] have proven the power of his method by applying it to the 
case of irrational parameters. Specifically, they have given the expansion of 

^J-4£,-i-5,-f-2e,i-35 
Table 3 shows their results, besides of ours. 



4 Final comments 

We have presented a direct and intuitive procedure to construct the er-expansion 
of any function, of one or several variables, of the hypergeometric class, when- 
ever the function can be written as a convergent power series of the variables. 
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Table 3 

Coefficients of the e-expansion of the 
Eg. 



generahzed hyper geometric function given in 



e order 


NumExp 


our method 


e° 


1.0000000000000000000 


1 


£^ 


-1.44555526747927 


-1.44555526747927 


£2 


3.9383879447 


3.938387944727 


£3 


-266.94735 


-266.9473544234 


8^ 


298.66 


298.666582673 






-56037.4042903 






-113001.082396 



The building blocks are the successive derivatives of Pochhammer and recip- 
rocal Pochhammer symbols particularized for the values of the parameters at 
e = 0. We have provided with closed explicit expressions for those derivatives 
to any order, although recurrence relations, given in Appendices A and B, can 
be equally used to evaluate them. Each coefficient of the e-expansion appears 
as a power series of the variable or variables that can be summed, to give a 
known function, only in very special cases. This fact, however, is not a serious 
drawback, in our opinion. In fact, the symbols commonly used to refer to most 
of the special functions are merely a shorthand to represent a series or an inte- 
gral. Nevertheless, what makes our procedure to be especially useful is its easy 
implementation to obtain numerical values of the coefficients. The examples 
given above show that our method allows one to calculate terms of higher order 
with a notable accuracy. Besides this, the proposed method is self-contained: 
it does not need auxiliary procedures like, for instance, differential reduction, 
even in the case of singular lower parameters. 

An important limitation to the use of the algorithm presented in this pa- 
per stems from the necessary convergence of the series representation of the 
function to be expanded in powers of e. This restricts the applicability of the 
method to certain domains of the variables. Obviously, one may have recourse 
to analytic continuation. For instance, in the very common case of q+iFq one 
could use, for |z| > 1, the relation [TBI, Eq. (10)] [21 Sec. 5.3, Eq. (3)] 



q+lFq 



'ax, 



/3. 



; Hq 



"h, {1 + tt/i - /3fc}fc=l,...,g 



r(ai)---r(ag+i; 
/ 

g+i-^ q 



Olh) 



j^{l + tt/i — 0.k\k=\,...,q+\\ki=h 



(50) 
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The discussion of the analytic continuation of all functions of the hypergeo- 
metric class is, however, out of the scope of this paper. 
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A Derivatives of the Pochhammer symbol (a)m with respect to its 
argument a 



We have introduced, in Eq. ([6]), the notation 



As (a)m is a polynomial of degree m in a, 

pW(c,) = o for k>m. (A.2) 



We assume in the rest of this appendix that m>k. A generating function for 
the P^^a) can be immediately obtained from [211 Sec. 6.2.1, Eq. (2)] 

oo 

E (a)m(-t)™/m! = iFo(a; ; -t) = (1 + t)-", |t| < 1. (A.3) 

m=0 



Derivation, k times, with respect to a gives 

oo 

E k\ {-tr/m\ = {-I f (1 + t)-- [ln(l + t)]' , |t| < 1. (A.4) 

m=0 

Then, we have 



/ 1 \ k—rn am 



t=0 



15 



k—m m 



E 

1=0 



Jm—l 



dt 



m—l 



[ln(l + t)f 



(A.5) 



i=0 



Obviously, 



Ql 



(A.6) 



i=0 



On the other hand, as [ln(l + t)]'^ is the generating function of the Stirhng 
numbers ^ Eq. 26.8.8], 



(A.7) 



[ln(l + t)]'' = k\J2 k) e/n\, \t\ < 1, 

n=k 



we get 



d 



m—l 



k\ s{m — l,k) . 



(A.8) 



i=0 



Substitution of Eqs. f lA.6P and (\A.8\i in ( 1A.5|) gives 

/ \ 

P«(«) = (-ir-'=5:(-iy(^ Js(m-/,fc)(a), for m > k. (A.9) 



For numerical computation, however, it may be preferable to make use of the 
recurrence relation 



= {a + m) V^J;\a) + V^-'\a), k > 0, 



(A.IO) 



with starting values 



(A.ll) 



B Derivatives of the reciprocal Pochhammer symbol with 
respect to its argument /3 



The notation 



1 1 



k\ dl3^ {/3)r 



(B.l) 
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has been proposed in Eq. ([6]) to represent the derivatives of the reciprocal 
Pochhamer symbol with respect to its variable. Very simple expressions for 
these can be easily obtained from the relation [221 Eq- 4.2.2.45] 



1 (-1 



/!(m-l-/)! /? + /' 



Direct derivation with respect to /3 in this equation gives 

m-l ( — lY 1 



provided /3 is different from a nonpositive integer, —N, such that < < m. 
To deal with this case (in the limit e — )• 0), we have defined, in Eq. (ITB]) . the 
regularized reciprocal Pochhammer symbol 

^ l3 = -N + be, 0<N<m. (B.4) 



whose derivatives have been denoted, in Eq. ( |T6l) . by 



From Eqs. ( IB. 21) and ( 1B.4P we obtain 



;eo/!(m-l-0! + / 



Repeated derivation with respect to (3 leads to 



For numerical implementation, one may use the recurrence relations 

QL'Ii(/3) = - QL';!^(/3)) + m), (B.8) 

with initial values 

Qr(/3) = 4,o, Q(^)(/3) = l/(/3U (B.9) 
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and, for the regularized reciprocal Pochhammer symbols when /3 = — + be, 
QtUd) = - + m), m>N, (B.IO) 

starting with 

Q^nIM = QP{P), = rn>N. (B.ll) 
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